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PREFACE 


The  latest  versions  of  the  atmospheric  general  circulation  models 
(GCMs)  incorporate  many  physical  processes  that  govern  the  larger 
scales  of  atmospheric  circulation,  although  there  is  still  need  to  im- 
prove their  fidelity.  This  need  has  been  considered  to  be  important 
with  respect  to  the  parameterization — or  physical  representation — of 
sub-grid  scale  turbulent  processes , and  specifically  the  parameteriza- 
tion of  the  fluxes  within  the  surface  boundary  layer. 

This  report  reviews  the  parameterization  techniques  currently  be- 
ing used  in  the  GCMs  and  suggests  a series  of  experiments  to  test  var- 
ious parameterization  schemes  in  the  Rand  two-level  GCM.  It  is  hoped 
that  the  findings  from  these  tests,  when  completed,  will  help  us  to  de- 
vise a scheme  that  represents  the  boundary  layer  processes  more  realis- 
tically. 

This  work  has  been  done  as  part  of  the  Rand  Climate  Program,  under 
support  from  the  Defense  Advanced  Research  Projects  Agency.  Related 
Rand  publications  on  the  dynamics  of  climate  are  R-877-ARPA,  A Doaimen- 
tation  of  the  M-intz-Avakawa  TixJO-leveZ  Atmospheric  General  Circulation 
Model  (December  1971)  ; R-1511-ARPA,  Pkmerical  Experiments  on  the  Com- 
putation of  Ground  Surface  Temperature  in  an  Atmospheric  Circulation 
Model  (May  1974);  and  R-1653-ARPA,  A Survey  of  the  Research  on  the 
Treatment  of  the  Planetary  Boundary  Layer  in  Atmospheric  Circulation 
Models  (in  preparation). 


Published  by  The  Rand  Corporation 


SUMMARY 


The  thermal  and  dynamic  interaction  between  the  atmosphere  and 
underlying  surface  occurs  through  the  boundary  layer.  While  these  in- 
teractions have  generally  been  ignored  for  short-term  prediction  of 
large-scale  atmospheric  circulations,  they  are  quite  important  for  long- 
range  forecasting  and  in  studies  related  to  general  circulation  and 
climate  dynamics.  In  this  case  a consideration  must  be  given  not  only 
to  the  supply  of  energy  but  also  to  the  dissipation  of  kinetic  energy, 
as  well  as  to  the  vertical  transport  of  heat  and  moisture  within  the 
boundary  layer.  Since  the  intensity  of  small-scale  processes  is  af- 
fected by  large-scale  processes,  incorporation  of  the  boundary-layer 
dynamics  constitutes  an  essential  part  in  studying  the  physical  prin- 
ciples of  long-range  forecasting  of  large-scale  processes  and  climatic 
changes . 

Numerical  general  circulation  models  have  been  adopted  as  impor- 
tant tools  for  understanding  the  physical  basis  of  climatic  changes. 
These  models  vary  greatly  in  many  respects,  especially  with  respect  to 
the  degree  of  detail  in  their  treatment  of  those  physical  processes 
that  cannot  be  resolved  by  the  grid  spacing  of  the  model  viz  the 
boundary-layer  processes. 

The  boundary-layer  parameterization  in  a general  circulation  model 
is  usually  dependent  on  the  vertical  resolution  of  the  model  and  is 
related  to  the  determination  of  four  factors  in  terms  of  the  variables 
predicted  by  the  model: 

(a)  surface  fluxes  of  momentum,  heat,  and  moisture, 

(b)  vertical  profiles  of  the  turbulent  fluxes  within  the  boundary 
layer, 

(c)  height  of  the  boundary  layer, 

(d)  vertical  velocity  at  the  top  of  the  boundary  layer. 

There  are  several  approaches  to  the  boundary-layer  parameteriza- 
tion aimed  at  incorporating  the  boundary— layer  processes  in  general 
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circulation  models.  Some  of  these  techniques  are  based  on  so-called 
K-theory  while  others  are  based  on  the  similarity  theory.  Considering 
the  varying  degrees  of  vertical  resolution  of  different  general  circu- 
lation models,  the  determination  of  the  surface  fluxes  is  perhaps  the 
most  important  aspect  of  the  boundary-layer  parameterization  in  a gen- 
eral circulation  model. 

At  present  there  is  not  sufficient  evidence  to  determine  which 
particular  boundary-layer  parameterization  scheme  is  most  satisfactory. 
Hence  it  is  suggested  that  a systematic  sensitivity  test  of  various 
schemes  be  carried  out  to  determine  the  best  approach. 
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TOTRODUCTION 


PARAMETERIZATION  AND  NEED  FOR  IT 

In  setting  up  numerical  models  of  atmospheric  processes,  the  scales 
of  primary  interest  are  included  within  the  available  resolutions  of 
grid  size.  For  example,  for  large-scale  motions  the  domain  of  interest 
is  the  entire  globe,  and  the  large-scale  processes  themselves  can  be 
resolved  by  prescribing  a horizontal  grid  width  of  a few  hundred  km  and 
a vertical  grid  mesh  size  of  50  m in  the  lower  atmosphere  and  several 
kilometers  in  the  free  atmosphere.  However,  resolutions  of  such  magni- 
tude cannot  account  for  significant  interactions  with  the  scales  that 
have  been  truncated.  Whereas  such  scales  have  been  neglected  for  short 
periods  a day  or  so)  , for  longer  periods  they  are  considered  to  be 
quite  important,  both  as  a source  as  well  as  a sink  of  energy  for  the 
explicitly  resolved  scales.  And  since  these  processes  (called  sub— grid 
or  sub-resolution  processes)  cannot  be  explicitly  handled  in  the  rele- 
vant numerical  model,  their  effects  are  usually  incorporated  by  ex- 
pressing their  statistical  effect  on  the  large-  (explicitly  resolved) 
scale  processes  in  terms  of  large-scale  parameters.  This  technique  of 
interactions  between  resolvable  and  unresolvable  processes  is 
called  f and  it  involves  empirically  related  parameters 

which  are  determined  either  on  the  basis  of  observations  or  from  theo- 
retical considerations. 

PROBLEMS  IN  PARAMETERIZATION 

Perhaps  the  most  important  problem  in  large-scale  modeling  is  to 
first  identify  which  sub-grid  scale  processes  are  of  importance.  This 
job  is  complicated  because  in  the  atmosphere  various  processes  (on 
different  scales)  tend  to  occur  simultaneously,  making  it  difficult  to 
isolate  a particular  scale  process  and  determine  its  role  meaningfully, 
especially  with  respect  to  its  interaction  with  the  resolved  scale 
processes.  In  view  of  the  complexities  and  interdependence  of  atmos- 
pheric processes  on  different  scales  and  also  the  lack  of  knowledge 
about  their  dynamics,  parameterization  of  some  sub-grid  scale  processes 
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is  generally  intuitive  until  the  scheme  is  actually  included  and  tested 
in  a large-scale  numerical  model.  Since  the  parameterization  techniques 
are  not  unique  and  can  be  approached  in  many  different  ways,  it  is  pos- 
sible to  get  involved  in  a lot  of  details;  this  should  be  avoided  un- 
less, of  course,  the  fundamental  physics  is  involved. 

The  first  and  most  critical  step  in  developing  a parameterization 
technique  is  to  be  able  to  define  the  problem  precisely.  It  is  impor- 
tant to  determine  which  of  the  several  different  smaller  scale  processes 
(present  simultaneously)  have  important  interactions  with  the  large- 
scale  process.  This  can  be  done  by  a combination  of  steps; 

o We  can  use  all  the  available  observations  of  the  relevant  sub- 
grid scale  process  to  obtain  a realistic  description  of  the  phenomena 
and  thus  determine  universal  values  for  empirical  constants. 

o We  can  define  parameterization  hypotheses  on  the  basis  of  ob- 
served physical  features.  These  hypotheses  can  be  tested  by  incorpor- 
ating the  parameterization  in  an  actual  general  circulation  model  (GCM) 
and  calibrating  it  until  a realistic  simulation  of  the  "modeled”  pro- 
cess is  attained.  For  this  purpose,  all  existing  explicit  models  of 
the  sub-grid  scale  processes  can  be  used  to  get  a better  insight  into 
the  physics  and  dynamics  of  the  process.  This  is  especially  true  in 
the  parameterization  of  the  atmospheric  boundary  layer  in  a GCM. 

We  have  seen  that  the  atmosphere,  which  physically  is  a continuum, 
is  treated  in  numerical  models  as  a fluid  with  finite  numbers  of  de- 
grees of  freedom.  The  number  of  degrees  of  freedom  is  determined  in 
the  horizontal  by  a distance  d between  discrete  grid  points  and  in  the 
vertical  by  a number  of  levels  n.  In  most  large-scale  models,  d ^ 

300  km  and  n ^ 20.  As  a consequence,  the  theoretical  model  equations 
cannot  describe  sub-grid  scale  boundary- layer  turbulent  processes,  and 
it  becomes  necessary  to  develop  certain  relations  for  stresses  and 
turbulent  heat  and  moisture  fluxes  in  terms  of  external  (to  the  bound- 
ary layer)  variables  resolved  and  computed  by  the  model. 

The  interaction  of  the  atmosphere  with  the  underlying  surface 
takes  place  via  the  atmospheric  boundary  layer  within  which  the  model 
equations  should  incorporate  the  predominantly  vertical  turbulent 
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fluxes  of  momentuin,  heat,  and  moisture.  These  fluxes  are  usually  large 
in  the  vicinity  of  the  underlying  surface  and  tend  to  decrease  with 
height.  Considering  the  vertical  resolution  of  current  global  circula- 

I'^odels , it  appears  desirable  to  neglect  the  detailed  vertical 
structure  of  the  boundary  layer  and  restrict  parameterization  of 
boundary-layer  processes  to  the  incorporation  of  near-surface  values 
for  horizontal  components  of  momentum,  heat,  and  moisture  fluxes  in  the 
model  equations.  Thus  the  parameterization  methods  should  be  able  to 
calculate  these  fluxes  in  terms  of  the  large-scale  output  of  the  model 
as  well  as  the  parameters  of  the  underlying  surface.  It  may  be  noted 
that  most  of  the  parameterization  techniques  used  in  GCMs  are  based  on 
the  local-  considerations  of  the  boundary  layer,  whereas  the  GCMs  them- 
selves deal  with  data  averaged  for  horizontal  distances  of  the  order 
of  several  hundreds  of  km.  While  the  effects  of  such  averaging  is  not 
likely  to  be  very  large  over  oceans  (Fleagle  et  al.,  1967),  they  can 
be  very  significant  over  land. 

Thus  the  purposes  of  this  report  may  be  stated  as: 

1.  To  discuss  the  interactions  between  the  atmospheric  boundary 
layer  and  the  large-scale  atmospheric  circulation. 

2.  To  discuss  the  parameterization  techniques  currently  being 
used  in  GCMs. 

3.  To  discuss  the  testing  of  various  parameterization  schemes  in 
a specific  GCM  in  order  to  estimate  the  advantages  and  deficiencies  of 
various  approaches. 
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II.  ATMOSPHERIC  BOUNDARY  LAYER  AND  LARGE-SCALE  ATMOSPHERIC  CIRCULATIONS 


It  is  generally  acknowledged  that  the  development  of  research  on 
atmospheric  boundary-layer  (b.l.)  dynamics  has  a direct  bearing  on  many 
problems  relating  to  micrometeorological  processes  within  the  b.l.  For 
example,  the  turbulent  state  of  the  boundary  layer  plays  a decisive  role 
in  the  phenomenon  of  the  diffusion  of  atmospheric  pollutants,  and  also 
the  dissipation  of  the  kinetic  energy  of  the  atmosphere  which  occurs  to 
a considerable  extent  in  the  boundary  layer.  It  is  through  the  boundary 
layer  that  there  is  thermal  and  dynamic  interaction  between  the  atmos- 
phere and  the  underlying  surface.  Thus  a knowledge  of  various  features 
of  the  b.l.  is  essential  for  numerical  weather  forecasting  and  also  for 
simulating  the  dynamics  of  climate.  Whereas  the  dissipative  role  of 
boundary-layer  turbulence  is  generally  ignored  for  the  shoi^t-range  fore- 
casting of  ta'pge-scaZe  atmospheric  circulations,  this  cannot  be  done 
in  the  case  of  long-range  forecasting,  general  circulation  studies,  and 
studies  of  climate.  This  is  because  in  the  latter  cases  (unlike  short- 
range  forecasting)  consideration  must  be  given  not  only  to  the  supply 
of  energy,  but  also  to  dissipation  of  kinetic  energy  within  the  b.l. 
and  to  the  vertical  transport  of  heat  and  moisture  through  the  boundary 
layer.  Since  the  intensity  of  micrometeorological  processes  is,  in 
turn,  affected  by  the  large-  (macro-)  scale  processes,  the  study  of 
b.l.  dynamics  must  constitute  an  essential  part  in  studying  the  phys- 
ical principles  of  long-range  forecasting  of  large-scale  phenomena 
and  the  simulation  of  climate. 

We  attempt  to  delineate  the  -interaction  between  the  b.l.  and  the 
free  large-scale  atmosphere  through  answers  to  some  specific  questions. 

Q.  What  is  the  basic  information  that  relates  the  large-scale 
processes  to  the  boundary r: layer  processes? 

The  interaction  between  the  atmosphere  and  the  underlying  surface 
occurs  through  the  b.l.,  and  the  vertical  transport  of  momentum,  heat, 
and  moisture  appears  to  be  the  basic  mechanism  of  exchanges  between 
the  atmosphere  and  the  earth.  Whereas  the  small-scale  turbulent  fluxes 
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of  the  quantities  are  continuous  and  reach  their,  maximum  xvithin  the 
b.l.,  they  are  intermittent  in  the  free  atmosphere,  beyond  the  b.l. 

The  exchange  between  the  b.l.  and  the  free  atmosphere  takes  place  mainly 
through  regular  vertical  transports  caused  by  evolution  of  meteorolog- 
ical fields  on  synoptic  scales  and  frictional  convergence  in  the  bound- 
ary layer.  Also,  the  exchange  caused  by  meso-scale  cumulus  convection, 
although  relatively  random,  occurs  through  the  boundary  layer  and  is 
essential  for  larger  scale  atmospheric  dynamics.  Figure  1 shows  the 
temperature  spectrum  in  the  surface  boundary  layer.  It  can  be  seen 
that  regions  of  the  spectrum  corresponding  to  the  small-scale  turbu- 
lence and  relatively  large  period  variation  of  temperature  are  divided 
by  a gap  of  periods  of  the  order  of  10  minutes  (Koleshnikova  and  Monin, 
1965).  In  other  words,  the  interaction  between  the  atmosphere  and  the 
underlying  surface  occurs  in  two  stages.  First  there  are  fluxes  through 
small-scale  turbulence  from  the  earth  to  the  boundary  layer,  and  at  the 
next  stage,  these  are  transformed  from  the  boundary  layer  to  the  free 
atmosphere  by  comparatively  regular  vertical  motions. 

Thus  the  basic  information  required  for  describing  the  large-scale 
processes  (in  terns  of  the  boundary-layer  processes)  relates  to  the 


Fig.  1 The  temperature  snectrum  in  the  surface  boundary  layer 
(Koleshnikova  and  Monin,  1965) 
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turbulent  fluxes  of  momentum,  heat,  and  moisture  near  the  underlying 
surface. 


Q.  Uhat  -is  the  important  link  between  the  b.l,  and  the  free  at- 
mosphere ? 

One  important  link  between  the  boundary  layer  and  the  free  atmos- 
phere is  the  vertical  velocity  at  the  top  of  the  boundary  layer.  Since 
boundary  layers  are  of  limited  vertical  extent,  it  is  vertical  velocity 
that  transports  momentum,  heat,  and  moisture  out  of  the  boundary  layer 
into  the  free  atmosphere.  This  velocity,  typically  of  the  order  of  mm/sec, 
is  largely  produced  by  frictionally  modified  flow  near  the  surface  and 
has  a significant  effect  on  the  production  of  cloud  and  the  general  de- 
velopment of  large-scale  circulations.  However,  the  mechanism  by  which 
the  boundary  layer  and  the  free  atmosphere  interact  is  not  understood 
adequately.  It  is  necessary  to  study  temporal  Variations  of  the  bound- 
ary layer  to  understand  the  mechanism  by  which  a previously  unstable 
boundary  layer  becomes  stable  or  vice  versa.  This  could  occur  by  ad- 
vection  and/or  by  entrainment.  The  latter  term  is  more  applicable  near 
the  top  of  the  boundary  layer  where  the  intensity  of  turbulence  and/or 
stability  of  the  air  above  the  boundary  layer  determines  the  extent  of 
incorporation  of  the  free  atmosphere  air  with  that  from  the  boundary 
layer  (or  vice  versa) . 

The  pressure  that  determines  the  gradient  wind,  and  the  precipi- 
tation which  falls  through  the  atmosphere  to  the  surface,  are  the  other 
links  between  the  free  atmosphere  and  the  b.l. 

Q.  Uhat  is  the  effect  of  underlying  surface  roughness  character- 
istics on  the  large-scale  atmospheric  circulations  I 

The  planetary  boundary  layer  is  defined  as  the  layer  of  air  in 
which  surface  frictional  effects  are  significant.  Studies  of  the  b.l. 
in  micrometeorology  have  examined  the  variations  in  momentum  flux  due 
to  spatial  variations  in  surface  roughness  characteristics  in  great 
detail.  Elliot  (1958) , who  studied  the  effect  of  abrupt  changes  in 
surface  roughness  on  momentum  flux,  suggested  that  an  internal  bound- 
ary layer  separates  the  flow  that  is  in  equilibrium  with  the  upwind 
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surface  from  the  flow  that  is  approaching  equilibrium  with  the  dovniwind 
surface,  Panofsky  and  Townsend  (1964)  generalized  Elliot's  model  and 
also  came  to  the  conclusion  that  there  is  a zone  of  transition  (in  the 
vertical)  that  separates  flow  above  that  has  not  "felt"  the  new  surface, 
from  flow  below  that  is  nearly  in  equilibrium  with  the  new  surface 
roughness.  In  addition,  various  other  studies  have  clearly  demonstrated 
the  importance  of  incorporating  variations  in  surface  roughness  in  local 
micro-meteorological  studies. 

There  is  considerable  interest  in  the  mechanisms  of  energy  trans- 
fer in  the  atmosphere  due  to  increased  interest  in  long-range  weather 
forecasting.  Lettau  (1959) , by  estimating  the  mean  reservoir  of  atmos- 
pheric mechanical  energy  and  the  global  mean  of  energy  dissipation,  has 
shown  that  there  is  a renewal  period  for  kinetic  energy  of  about  three 
days.  This  would  suggest  that  any  long-term  prediction  scheme  for  a 
period  greater  than  one  or  two  days  must  include  a proper  specification 
of  mechanisms  for  frictional  dissipation.  As  a consequence  it  is  nec- 
essary to  obtain  a fairly  precise  knowledge  of  both  spatial  and  temporal 
variations  of  dissipation  effects  on  (large)  scales,  Lettau  has  at- 
tempted, both  empirically  (1959)  and  theoretically  (1962),  to  solve  the 
general  problem  of  including  surface  fluxes  in  the  large-scale  systems 
by  relating  them  to  the  external  parameters  controlling  the  flow.  He 
introduced  the  concepts  of  geostrophic  drag  coefficient  and  surface 
Rossby  number,  both  of  which  depend  on,  among  other  factors,  the  sur- 
face roughness  length.  Kung  (1963)  subsequently  used  Lettau 's  model 
to  estimate  climatological  patterns  of  energy  dissipation  over  the  en- 
tire Northern  Hemisphere.  His  calculations  were  based  on  his  assess- 
ments of  the  surface  roughness  parameter,  which  in  turn  were  based  on 
land  use  statistics  at  various  (360)  locations.  Holopainen  (1963)  has 
also  estimated  frictional  dissipation  over  the  British  Isles  for  a 3- 
month  period. 

On  a global  scale  the  surface  characteristics  can  vary  over  a wide 
range,  as  over  mountains  of  varying  sizes  and  shapes,  forests,  etc. 

As  indicated  by  Sawyer  (1959) , large-scale  circulations  of  the  atmos- 
phere are  affected  by  the  characteristics  of  the  underlying  surface 
through  different  mechanisms  which  depend  on  the  physical  extent  of 


surface  characteristics.  For  example,  large-scale  motion  can  be  affected 
by: 

o Surface  roughness  over  land  features  whose  dimensions  may  range 
from  a few  cm  to  a few  hundred  meters,  as  well  as  those  rang- 
ing from  a few  km  to  several  hundreds  of  km. 

o Differences  between  the  roughness  characteristics  of  land  and 
ocean. 

In  each  of  these  cases  it  may  be  necessary  to  consider  the  effect  of 
local  surface  characteristics.  For  example,  in  a global  circulation 
study,  it  may  be  desirable  to  incorporate  the  effect  of  a mountain  clb 
welt  as  the  local  "roughness"  characteristics  of  the  mountain  itself. 

Q.  How  does  the  surface  drag  within  the  boundary  layer  modify 
the  large-scale  atmospheric  circulations!. 

Charney  and  Eliassen  (1949) , who  were  the  first  to  pose  this  ques- 
tion, have  described  the  principal  dynamical  mechanism  involved.  It 
is  argued  that  motion  with  the  friction  (surface)  layer  occurs  with 
approximate  balance  among  the  horizontal  pressure  gradient , Coriolis 
acceleration,  and  the  eddy  stress.  This  results  in  the  Ekman  drift 
across  the  isobars  toward  low  pressure,  which  causes  horizontal  conver- 
gence of  air  within  the  friction  layer  into  regions  of  cyclonic  vortic— 
ity  and  divergence  from  regions  of  anticyclonic  vorticity.  Conse- 
quently, a compensating  net  divergence  must  occur  in  the  free  atmo-' 
sphere.  There  is  a positive  correlation  between  divergence  and  vortic- 
ity in  the  free  atmosphere  which  leads  to  deceleration  of  the  winds. 

This  theory,  of  course,  would  fail  in  case  of  deep  convection  where 
the  effective  eddy  stress  is  not  confined  to  a relatively  shallow  layer, 
and  also  over  the  oceans  where  the  eddy  stress  in  the  free  atmosphere 
may  be  comparable  with  the  surface  stress  (Sheppard  et  al.,  1952). 

Over  land  the  theory  is  quite  useful  since  eddy  stress  in  the  surface 
layer  is  much  greater  than  that  in  free  atmosphere.  This  theory  has 
been  used  to  determine  the  most  important  link  between  the  b.l.  and 
free  atmosphere,  namely  the  vertical  velocity  (W^^)  at  the  top  of  the 
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b.l,  in  terms  of  free  atmospheric  variables,  such  as  geostrophic  vor- 
ticity  at  the  top  of  the  b.l.  Mahrt  (1974)  has  extended  Charney  and 
Eliassen  s theory  to  depict  as  a function  of  geostrophic  relative 
vortlcity,  drag  coefficient,  and  boundary-layer  depth. 


Q.  To  what  ho'ight  does  the  b.l.  extend  vnto  the  atmosphere  ccnd 
what  aj'e  the  factors  which  determine  this  heightl 

The  b.l.  thickness  varies  not  only  in  time  but  also  in  space.  The 
study  of  boundary- layer  thickness  and  its  variation  from  hour  to  hour, 
from  day  to  day,  and  from  season  to  season  has  become  the  subject  of 
extensive  research.  The  boundary-layer  thickness  depends  on  many  fac- 
tors, including  the  following; 

o stability  of  the  surface  layer  which  in  turn  depends  upon  the 
time  of  the  day, 

o entrainment  of  free  atmosphere  air  into  the  boundary  layer 
near  its  top, 

o large-scale  vertical  motions  at  the  top  of  the  boundary  layer, 

o horizontal  advective  changes  near  the  boundary-layer  top  (which 
control  development  of  boundary  layer  over  sea) , 

o radiative  heat  fluxes  and  latent  heat  effects  which  are  sig- 
nificant when  there  are  clouds  present  within  the  boundary 
layer, 

o the  lapse  rate  of  potential  temperature  in  the  air  above  the 
boundary  layer  and  the  intensity  of  the  inversion  in  case  of 
mixed  layers . 


Figure  2 shows  the  diurnal  variation  of  the  mean  boundary-layer 
thickness  under  strongly  convective  daytime  conditions.  It  can  be  seen 
that  it  varies  from  a few  tens  of  meters  at  night  to  nearly  2000  m dur- 
ing the  afternoon.  By  late  afternoon  there  is  no  further  thickening 
of  the  boundary  layer,  and  around  sunset  the  boundary  layer  becomes 
quite  shallow  and  grows  slowly  through  the  night  to  complete  the  diur- 
nal cycle.  Thus  there  are  significant  discontinuities  in  the  b.l. 
depth  near  sunrise  and  sunset.  A detailed  quantitative  treatment  of 
determining  the  height  of  the  b.l.  has  been  discussed  by  Bhumralkar 
(1974)  in  a separate  survey  report. 
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Fig.  2 — The  mean  boundary-layer  thickness,  <h(t)>,  deduced 
for  the  O’Neill  data  and  plotted  with  standard  errors  as 
functions  of  time  of  day,  t,  in  Mean  Solar  Time 
(Carson,  1973) 


From  the  above  it  is  clear  that  treatment  of  the  boundary  layer 
for  studies  of  the  general  circulation  of  the  atmosphere  presents  much 
more  formidable  problems  than  those  encountered  in  smaller  (micro/meso/ 
synoptic)  scale  studies.  The  most  important  problems  pertain  to  the 
wide  range  of  variations  in  the  b.l.  characteristics  over  the  entire 
globe  and  the  lack  of  a "universal”  theory  to  study  these. 
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III.  GOVERNING  EQUATIONS  OF  THE  GENERAL  CIRCULATION  MODEL 
RELEVANT  TO  THE  BOUNDARY-LAYER  PROCESSES 


INTRODUCTION 

The  dynamical  (numerical)  general  circulation  models  (GCMs)  have 
been  developed  with  the  immediate  objective  of  understanding  the  gen- 
eral circulation  of  the  atmosphere.  They  have  also  been  adopted  as 
important  tools  for  the  longer  range  objective  of  understanding  the 
physical  basis  of  climatic  changes. 

During  the  last  decade  or  so  there  has  been  significant  progress 
in  the  physical  and  theoretical  basis  for  the  development  of  numerical 
modeling  of  the  general  circulation  due  to  our  increased  understanding 
of  physical  processes  in  the  atmosphere.  However,  most  of  these  pro- 
cesses have  been  incorporated  in  GCMs  through  parameterization  tech- 
niques only.  The  type  and  degree  of  parameterization,  of  course,  de- 
pends on  the  horizontal  and  vertical  resolution  of  the  particular  model. 
For  example,  on  one  hand  the  most  simple  models  are  the  vertically  and 
zonally  averaged  models  (Sellers,  1973;  Budyko,  1969)  which  parameter- 
ize the  transport  of  heat  in  terms  of  mean  zonal  variables;  on  the 
other  hand,  there  are  models  of  very  high  resolution  which  are  capable 
of  resolving  the  details  of  cyclone-scale  motions.  However,  even  the 
models  of  high  resolution  are  required  to  parameterize  certain  physical 
processes  through  empirical  or  statistical  representation  of  nonresolved 
(sub-grid  scale)  processes  in  terms  of  the  resolved  parameters. 

All  GCMs  are  based  on  the  fundamental  hydrodynamical  equations 
that  govern  the  large-scale  behavior  of  the  atmosphere.  The  governing 
equations  of  a GCM  consist  of  the  momentum,  hydrostatic,  thermodynamic 
energy,  continuity  of  mass,  and  water  vapor  continuity  equations.  They 
are : 


3V  ^ 

— + V-VV  + 
dt 


9z 


+ 2J^xV+-Ap  = F 


(1) 


iB.  + Pg  = 0 , 
3z 


(2) 
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|i  + ?.ve  + w If  = 


(3) 


9p  3 

~ + v.pV  + — (wp)  = 0 . 


(4) 


If  If  = 

P = PRT  . (6) 

Here  F,  , and  represent,  respectively,  the  sources  and  sinks 
of  momentum,  heat,  and  moisture  due  to  various  physical  processes  in 
the  atmosphere.  The  system  of  governing  equations  (1  ^ 6)  is  closed 
by  parameterizing  F,  , and  S in  terms  of  the  large-scale  Cdependent) 
variables  V,  9,  q,  and  p,  GARP  Publications  Series  No.  14  (June  1974) 
has  compiled  the  current  status  of  global  models  and  has  given  details 
of  the  GCMs  that  are  in  use  all  over  the  world.  It  is  evident  that  the 
models  vary  considerably  in  many  respects,  particularly  with  regard  to 
the  degree  of  detail  in  their  treatment  of  the  sub-grid  scale  physical 
processes.  At  present  there  is  neither  sufficient  experimental  nor 
theoretical  evidence  to  determine  which  particular  treatment  (param- 
eterization) of  the  unresolved  physical  processes  included  in  F,  S , 
and  is  the  most  satisfactory. 

GOVERNING  EQUATIONS  RELEVANT  TO  THE  BOUNDARY  LAYER 

As  stated  earlier  (Sec.  I)  the  most  important  problem  in  large- 
scale  modeling  is  to  first  identify  which  sub-grid  scale  processes  are 
of  importance.  Here  we  attempt  to  identify  the  contribution  of  the 
b.l.  processes  to  the  large-scale  general  circulation  and  write  the 
governing  equations  of  the  GCM  as  if  other  physical  processes  were 
nonexistent . 

Physically,  the  term  F (Eq.  1)  represents  the  net  frictional 

force  consisting  of  the  frictional  drag  at  the  earth’s  surface  (F  ) 

->  s 

as  well  as  the  internal  friction  in  the  free  atmosphere  (F^)  . The 
former  is  the  part  related  to  the  boundary  layer.  The  flowchart  in 
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Fig.  3 shows  schematically  the  various  components  of  F (Box  numbers  3 

s 

to  5)  which  are  to  be  parameterized  in  a GCM  for  incorporating  the 
physical  effects  of  the  boundary-layer  friction. 


1 


Fig.  3 — Flow  diagram  showing  the  components  of  F , where 

s 

and  CF  ) represent,  respectively,  the  horizontal  and 
® V 

vertical  components  of  the  frictional  drag  at 
the  earth’s  surface. 


The  term  (Eq.  3)  represents  the  net  diabatic  heating  rate  and 
consists  of  three  parts: 


o S^;  heating/cooling  due  to  radiation, 

o ' l^^^ting  due  to  release  of  latent  heat  during  condensation 

o : heating/cooling  due  to  turbulent  heat  transports. 


It  is  evident  that  the  last  part  represents  boundary-layer  effects 

The  flow  diagram  (Fig.  4)  shows  the  components  of  S (box  numbers  5 to 

6 

9)  which  are  relevant  to  the  boundary-layer  processes  and  should  be 
parameterized  in  a GCM. 

The  terms  , and  are  evaluated  only  at  levels  within 

the  b.l.,  and  the  term  is  evaluated  at  all  levels  of  the  GCM. 

Finally,  the  term  (Eq.  5),  which  represents  the  net  moisture 
addition  rate,  consists  of  the  difference  between  the  evaporation  rate 
and  the  condensation  rate.  The  evaporation  can  take  place  both  from 
the  underlying  surface  as  well  as  clouds  and  precipitation  in  the  free 
atmosphere.  For  the  boundary-layer  processes,  we  will  be  mainly 
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Fxg.  4 -Flow  diagram  showing  components  of  S_,  where  S = radiation  effects 
(short  wave  + long  wave) ; S^=  consensation  in  free  atmosphere 
(cumulus  scale  4-  large  scale)  ; interaction  between  the 

b.l.  and  cumulus  scale  processes;  S^=  turbulent  transport  of 
heat  (boundary  layer  effects) ; convergence  of  vertical 

flux  of  sensible  heat;  horizontal  diffusion  of  sensible 

heat;  condensation  effect  when  fog/stratus  clouds  are 

present  within  the  boundary  layer;  3^^^=  radiation  cooling 
within  the  boundary  layer  when  it  contains  fog/stratus  clouds. 


concerned  with  the  evaporation  from  the  underlying  surface  as  well  as 
from  fog/stratus  clouds  which  may  be  present  within  the  boundary  layer. 
This  can  be  brought  about  by  the  turbulent  flux  of  moisture,  both  in 
the  vertical  and  horizontal  directions.  The  flow  diagram  (Fig.  5)  il- 
lustrates the  parts  of  whi-oh  ave  restated  to  the  boundary  layer ^ 

These  are  indicated  by  box  numbers  7 to  10. 

In  view  of  the  above  considerations,  we  can  rewrite  the  governing 
equations  for  momentum,  thermodynamic  energy,  and  water  vapor  for  the 
part  related  to  the  boundary  layer  as  follox^fs: 


dt 


+ (F  ) + (F  ) 

V H 


(7) 


■ ®DV  =DH  Sl  " ^DR  • 


dt 


(8) 
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Fig-  5 — Flow  diaeram  showi„ncr  the  co-^ronents  oh  , vrhere  P=  precipitation; 

cumulus  scale  precipitation;  P^=  large-scale  precipitation; 

E=  evaporation;  Eg=  boundary  layer  effects;  E^=  evaporation 

from  clouds  and  precipitation  in  free  atmosphere;  E^,=  vertical 

BV 

water  vapor  flux  through  the  b.l.  (evaluated  in  the  b.l.  only); 

E = horizontal  diffusion  of  water  vapor  x^ithln  the  b.l.  (evaluated 
Btl 

at  all  levels  in  the  GCM) ; E„^=  evaporation  from  fog/stratus 

Br 

clouds  within  the  boundary  layer.  (Note:  Interaction  between 

the  b.l.  and  cumulus  scale  processes  and  the  effect  of  precipitation 
falling  through  the  b.l.  should  also  be  considered.) 


dt 


"BV 


+ E 


BH 


- E 


BF 


(9) 


where  the  terms  on  the  right-hand  sides  of  Eqs.  (7) -(9)  have  been  iden- 
tified above  in  the  flow  diagrams  for  F,  Sq,  and  S . It  may  be  men- 
tioned here  that  in  these  equations  we  have  not  included  contributions 
due  to  interactions  between  the  b.l.  and  the  cumulus  scale  processes, 
and  consequently  the  terms  on  the  right-hand  side  are  strictly  for  the 
part  related  to  the  b.l.  processes  only. 

The  problem  of  parameterization  of  the  b.l.  processes  in  a GCM 
thus  reduces  to  the  determination  of  the  terms  F^,  F^, 

S E . E , and  E in  terms  of  the  large-scale  variables  provided 
DR’  BH’  BV’  BF 

by  the  GCM. 
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IV.  PARAMETERIZATION  OF  THE  BOUNDARY-LAYER  EFFECTS  IN  CURRENT 
ATMOSPHERIC  GENERAL  CIRCULATION  MODELS 


INTRODUCTION 

There  are  various  atmospheric  models  with  primitive  equations  that 
are  either  hemispheric  or  global  and  with  varying  degrees  of  vertical 
resolution.  The  ways  of  treating  the  various  physical  processes  differ 
considerably  depending  on  the  purpose  for  which  the  models  are  designed. 
Most  present  models  are  applied  to  weather  forecasting  for  1 to  5 days, 
while  some  are  being  developed  as  tools  to  study  atmospheric  changes 
beyond  a few  weeks  or  longer,  i.e.,  to  study  climatic  variations.  In 
this  report  we  consider  this  latter  type  of  GCMs  only.  These  are  the 
GCMs  developed  by 

o Geophysical  Fluid  Dynamics  Laboratory  (GFDL) , NOAA,  USA, 
o National  Center  for  Atmospheric  Research  (NCAR) , USA, 
o The  Rand  Corporation,  USA, 

o Goddard  Institute  for  Space  Studies  (GISS) , NASA,  USA, 
o University  of  California  at  Los  Angeles  (UCLA),  USA, 
o Meteorological  Office  (UK-I)  and  Universities  Modeling  Group 
(UK-II),  UK, 

o Hydrometeorological  Research  Center,  USSR, 

We  will  now  consider  specific  formulations  which  have  been  used 
in  various  GCMs  to  determine  the  terms  on  the  right-hand  side  of  Eqs. 
(7)-(9).  Whereas  all  the  GCMs  considered  in  this  report  evaluate  the 
vertical  eddy  fluxes  of  momentum,  heat,  and  moisture  (represented  by 
Fy>  ^BV’  tespectively)  , only  some  incorporate  horizontal  sub- 
grid scale  mixing  (F^^,  and  through  use  of  the  nonlinear  eddy 

diffusion  coefficient  proposed  by  Smagorinsky  (1963).  The  importance 
of  {.nteTaotion  between  the  b.l.  and  cumulus  scale  processes  has  been 
emphasized  both  observationally  (Gray,  1972)  and  theoretically  (Ara- 
kawa  and  Schubert,  1974).  However,  this  Interaction  (term  Fig-  4) 

is  a complicated  problem  requiring  much  further  research,  and  we  will 
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not  consider  it  in  our  report.  The  other  physical  effects  of  b.l.  which 
have  not  been  generally  included  in  GCMs  pertain  to  condensational  and 
radiative  effects  (terms  and  Eq.  8)  within  the  b.l.  Again 

these  have  been  included  only  in  the  UCLA  GCM,  which  considers  boundary 
layers  with  stratus  clouds  as  well  as  fog. 

From  the  above  discussion  it  appears  that  the  terms  that  have  been 
parameterized  to  incorporate  b.l.  effects  in  GCMs  are  those  which  repre- 
sent the  vertical  turbulent  fluxes  of  momentum  (F^) , heat  , and 

moisture  (E^^) , which  of  course  include  surface  fluxes.  Thus  each  of 
the  terms  F^,  and  E^^  consists  of  two  parts:  (1)  the  flux  (F  ) 

within  the  surface  layer  (usually  considered  to  be  a constant  flux 
layer)  extending  from  the  underlying  surface  to  a few  tens  of  meters 
above,  and  (2)  the  flux  (F^)  within  the  so-called  Ekman  layer  that  lies 
from  the  top  of  the  surface  layer  to  the  top  of  the  b.l.  Thus  we  have 


= F + F 
V o u 


Snn  = S + S 
DV  o u 


E = E + E 
BV  o u 


Any  b.l.  parameterization  in  a GCM  must  be  able  to  treat  the  b.l. 
processes  in  a physically  realistic  manner  and  relate  the  calculated 
variables  (within  the  GCM)  to  the  following  four  factors: 

(1)  the  surface  values  of  the  turbulent  fluxes  of  momentum,  heat, 
and  moisture; 

(2)  vertical  profiles  of  the  turbulent  fluxes  (as  well  as  of  the 
meteorological  parameters)  within  the  b.l.; 

(3)  the  height  of  the  b.l.; 

(4)  the  vertical  velocity  at  the  top  of  the  b.l. 


We  now  consider  the  determination  of  the  above  factors  by  the 
current  GCMs.  Considering  the  varying  degrees  of  vertical  resolution 


r^X8^ 

of  different  GCMs  (which  may  or  may  not  enable  determination  of  actual 
profiles  of  fluxes) , determination  of  the  surface  fluxes  is  perhaps  the 
most  important  aspect  of  b.l.  parameterization  in  GCMs. 

DETERMINATION  OF  SURFACE  FLUXES 

The  actual  method  of  parameterization  mainly  depends  on  whether  or 
not  the  b.l.  is  included  explicitly  in  the  numerical  model.  It  also 
depends  on  the  closeness  of  the  first  internal  GCM  level,  at  which  the 
complete  set  of  model-governing  equations  is  solved,  to  the  underlying 
surface.  For  example,  if  the  first  GCM  level  is  such  that  the  whole 
b.l.  is  placed  below  it,  then  only  the  surf ace/ near-surf ace  values  of 
the  turbulent  fluxes  of  momentum,  heat,  and  moisture  are  relevant  for 
the  b.l.  parameterization,  and  these  are  determined  diagnostically. 

This  can  be  achieved  by  the  use  of  either  some  resistance,  heat  trans- 
fer, and  moisture  transfer  laws  or  by  some  bulk  transfer  formulas. 

These  techniques  are  also  applicable  when  the  lowest  GCM  level  is 
placed  within  the  surface  layer  (or  constant  stress  layer)  at  a height 
of  about  ten  or  a few  tens  of  meters. 

Case  I.  Bulk  Type  Formulas 

The  surface  fluxes  are  calculated  by  the  bulk  formulas: 


— >■  1 I 

F = T = p C^V  V 
o o s D s ' s ' 


1 


S 

o 


= - p C C„  V 
s p H'  s 


(T. 


(11) 


E = E 
o o 


P (GW)C^ 
s E 


(q. 


(The  above  equations  are  solved  together  with  the  heat  budget  equation.) 

Here  t , Q , and  E are  the  surface  stress,  heat  flux,  and  moisture 
o o o ^ 

flux,  respectively.  V^,  T^,  q^ , and  are  wind  vector,  temperature, 
moisture,  and  density  at  the  level  s.  T^ , q^  are  temperature  and  mois- 
ture at  the  underlying  surface,  and  GW  is  the  ground  wetness  (0  1). 
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Case  II,  Mon in- Obukhov  Similarity  Theory  Formulas 

According  to  this  theory  the  surface  values  of  the  turbulent  fluxes 

are  given  by  a,  which  is  the  angle  between  the  surface  stress  (t^  ) and 

V , and  by  the  values  of  U . , Q , and  E . The  formulations  are 
s * ’ o 


c„|v 

s p H ' s 


A6 


(12) 


= - P (GW)C„|V  |Aq  . 
o s E s 

Here  A6 , and  Aq  are  the  differences  between  potential  temperature  and 

the  mixing  ratio,  respectively,  at  the  boundari-es  of  the  surface  layer. 

In  order  to  evaluate  t , Q , and  E in  GCMs  (from  Eqs.  (11)  and 

o o o 

(12)) the  knowledge  of  the  following  parameters  is  essential: 

o height  of  the  top  of  surface  layer  (for  Eq.  (11))  and  height 
of  the  b.l.  (for  Eq.  (12)); 

o ^s’  ^s’  *^s’  wind,  temperature,  and  moisture  at 

level  s ; 

o the  temperature  and  moisture  at  the  underlying  sur- 

face ; 

o GW,  the  wetness  measure  of  the  ground  (0  1)  ; 

o values  of  the  coefficients  C , C„,  and  C ; 

I)  n L 

o Z^,  the  roughness  parameter  (Eq.  (12)  only); 
o a (Eq.  (12)  only) , the  angle  between  the  surface  stress  and 


o (Eq.  (12)  only) , friction  velocity. 


Table  1 gives  a description  of  the  values  of  the  parameters  in- 
volved in  formulations  of  both  Case  I and  Case  II  and  used  in  various 
GCMs.  It  can  be  seen  that  out  of  eight  GCMs  considered,  four  use  bulk 
aerod3mamic  formulations  (Eq.  (11))  to  evaluate  surface  fluxes  whereas 
three  others  use  Eq.  (12) . The  technique  used  by  the  UCLA  GCM  is 


Table  1 


1.  CFDf.  (197J> 


Holloway  Si  Kanabe 
(1971) 


2.  NCAR  (1971) 


Kasahara  A 
Washington  (1971) 


RAND  (1971) 


Ciates  et  a.l . 
(1971) 


4,  GISS  (1974) 


Somerville  et  al, 
(1974) 


5.  UCLA  (1974) 


Arakawa  and  Mint^ 
0974) 


COMPU'CATION  OF  SURFACF  I'LUXKS  IN  VARIOUS  CCMS 


Fluxes 
Calculated 
from  Eq. 
01)  or 
iiq.  (12) 


Cnso  I 
T.q-  (11) 


Case  I 

K.(.  (U) 


No . of 
Levels 
Coor- 


Location  of 
Leve  1.  s 


Lowest  GCK 
Level /Layer 


Values  for  C /C  /(! 


1)  h'  li 


Case  I 

Eq.  (11) 


Case  I 

Eq.  (11) 


No  explicit 
level  s in 
the  GCM 

Lowest  cell 
level  = 3 km 


2/c  I No  explicit 
level  s in 
the  GCM 

Lowest  GCM 
level  = 800  mb 
- a « .75 
layer  = 600  mb 
= (’  = .50 


S)  ^ S.  ' 

,^-3 


(i)  2 X 10  (both  land  and  sea) 

(ii)  (C^)ocean  = 1.1  x lo’^ 

(i;^)  land  - 4.3  X lO"^ 


Cp  = Cj,  = 3 X 10-3 

Cj.  = 0.7  C to  reduce  evaporation 


No  explicit 
level  s in 
the  GCM 

Lowest  GCM 
level  ~ I km 
layer  ~ a = . 9 


For  neutral  .surface  Layer 

CC^)ocean  = min[(i  + .07 Iv  |)  I0" 
.0025] 

(C  )land  = .002  + .006  ( — — ) 
^ \5000  m/ 

also  C,  = C„ 


D 


H 


For  neutral  surface  Javer  .same  as  in 
RAND  GCM 

For  non-noutcal  condition.s  (for  all 
Rurface.s) 


f) 


(C  )aeuCral  


I. 


1-7  AT/|V^|' 
for  AT  c 0 (stable) 

“"D  " «’p)neutraJ  (1  + 

for  AT  ■>  0 (unstable) 

AT  = T,  - r 
b a 

also  Cr,  “ C,, 


3/o  I No  explicit 
level  s in 
the  GCM 

Lowe.st  GCM 
level  (1  = 8/9 
Layer  a = 7/9 


H 


V. 


and  have  a functional  dependence 

on  bulk  Richardson  number  and  (h/Z  ) 

, , o 

n IS  tlie  depth  of  b.  !.• 

(.45  m(l,ind 
is  roughness: I ^ 

^2.5  X 10  m(oceans) 


Determination  of 


All  three  calculated  from  prognostic 
equations  at  the  level  s 


f All  three  calculated  diagnostically  by 
equating  surface  fluxes  to  the  fluxes  In 
layer  from  s to  the  lowest  GCM  level  (at 
3 km)  (see  text) 


Sea:  speci- 
fied 

Land:  heat 
budget 
equation  at 
the  Surface 


Saturated  q at 


Same  way  as 
GFDI. 


obtained  by 
linear  extrap- 
olation of 
predic  ted 
velocities  at 
levels  1 & 3 
to  the  lower 
boundary 


•T^,  9g  arc  obtained  by 

equating  surface  fluxes 
of  sensible  heat  and 
moisture  to  the  same 
fluxes  from  s to  the 
lowest  GCM  level  (at 
o * .73 


Prognostic 

equatic^n 


Vg  obtained  by 

extrapolating 
ve.locities  in 
layers  8 .and  9 
to  the  lower 
bounda  ry 


Tg»  q^  same  as  in  RAND 
except  the  surface  fluxes 
are  equated  to  fluxes 
from  5 to  a = .9,  the 
lowest  GCM  level 


Prognostic 

equation 


GIDL  (1974):  18  levels  in 
the  vertical;  level  s at 
36  tn 


Same  as  CFDL 


Saturated  q at 


Saturated  q at 


1 

ro 

0 

1 


Table  1 (continued) 

COMPUTATION  OF  SURFACE  FLUXIiS  IN  VARIOUS  GCMS 


Model 

Fluxes 
Calculated 
f roTTi  Eq . 
(11)  or 

Eq.  (12) 

No.  of 
Levels 
Coor- 
dinate 
System 

Location  of 
Level  s 

Lowest  GCM 
Lcvel/Laycr 

Values  lor  C /C  /C 
D H E 

Determination  of 

V 

s 

L 

4 

% 

iemark-s 

6.  UK  I (19T2) 

(Similarity  theory) 

Corby  et  al . 

(1972) 

Case  II 
Eq.  (12) 

5/o 

No  explicit 
level  s in 
the  GCM 

Lowest  GCM 
level  at  1 kjji 

(also  the  top 
of  the  b-1.) 

"d  * S ' S 

0,3  X 10  ^ stable  n/6 

Land; 

[4.0  X lO"'^  unstable  r/lB 
To.  2 X 10  ^ stable  n/9 

Ocean: 

[2.0  X 10  unstable  0 

(1  is  the  angle  between  the  surface 
stress  and  the  wind  at  level  s.  ' 

roughness  parameter 
^0  ~ land;  1 cm  for  sea 

Computed  prog- 
nosticallv  at 
1 km. 

''«e  - (^4  - 

(“e>b 

Aq  = (q)^  - 

Uses  equivalent  potential 
temperature  (see  text) 

e 

where  9^  = equivalent  po- 
tential temperature. 

^ ^ ] 
• Aq 

Prognostic 

equation 

v75%  over, 
/ ocean  \ 

RH  4 

\ 40%  over / 
Tand  ' 

Formulation  for  0 Is 
0 

modified  in  order  to 
establish  a near  moist 
adiabatic  lapse  in 
lowest  km  over  typical 
oceans  (as  is  observed) 
and  at  the  same  time 
allow  for  a steeper 
lapse  rate  to  exist 
over  treated  land 
areas  (see  text). 

7,  UK  II  (1974) 

(Similarity  theory) 

(Only  ocean-covered 
planet  considered) 

Pearce  (1974) 

Case  II 
Eq.  (12) 

5/a 

No  explicit 
level  s in 
the  GCM 

I.owest  level 
at  850  meters 

■ ([Tf)’  ■ 

are  obtained  from  nomograms  (Clark. 
1970)  which  ace  plots  of  U^ , angle  a, 
C , C_  as  functions  S and  R , 

M.  h p 

Given  by  prog- 
nostic equa- 
tion at  850  m. 
(This  is 
required  to 
obtain  Q and 

V-  ° 

(i)  Surface 
stress  given 
by 

t “ 
o * 

(ii)  T , Q 
0 ^0 

and  E are 
zero  when 

|v^|  < 5 ~ 

T^»  are  referred  to 

the  lowest  GCM  level, 
i.e.,  at  850  m to  obtain 
|AT  = 

liq  = (q)^  - 

Prescribed 
as  function! 
of  latitude 

Saturated  q at 

S — stability  parameter 

= (g  ^ + eg4q)/.ff  [ |V^| 

AT  = T - T 
0 h 

R = surface  Rossby  No.  = 

l?J 
■ Ifp; 

0 = potential  temperature 
is  used  instead  of 
Temperature 

f corresponds  to  21'* 
latitude  between  21® 
and  equator 

8.  USSR  (1974) 

Fux-Rabinovich 

(1974) 

Case  11 
Eq.  (12) 

5/r. 

No  explic.U 
level  s in 
the  GCM 

Upper  limit 
of  the  b.l. 
= 1 . 5 km 
fixed 

S’  ''ll'  from 

nomograms  (same  as  in  UK  11  above 

(7)). 

Details  not  available 

4,  coordinate  in  ver- 
tical: takes  into 

account  the  global 
Orography 

g - hY„ 

T j V 1 

h:  b.l.  top  (first  GCM 

level) 

Y = lapse  rate  at  h - 
6®C/km 

22- 


discussed  below.  While  most  GCMs  use  Eq.  (12)  to  determine  the  surface 
fluxes,  UK-I  (Corby  et  al. , 1972)  has  modified  it  as  follows: 


1.  Since  UK-I  GCM  does  not  predict  surface  wind,  the  surface 
stress  is  prescribed  to  act  in  a different  direction  from  that  of  the 
wind  at  the  lowest  GCM  level  (i.e.,  lowest  free  atmosphere  wind).  As 
a consequence,  if  a is  the  angle  by  which  surface  wind  is  backed  from 
lowest  free  wind,  then  the  A and  c|)  components  of  surface  stress  are 
given  by  (instead  of  of  Eq.  (12)) 


T^o  “ cos  a - V sin  a) 

^(j)0  ""  ^sSl^s^  cos  a + U sin  a) 


(13) 


Also,  Cj^  and  a are  functions  of  stability  and  surface  roughness  with 
the  values  given  in  Table  1. 

2.  (of  Eq.  (12))  is  replaced  by 


- pc  C„ 

p H 


|(A0  + e^) 


(14) 


where  A = 0 for  stable  conditions.  This  result  is  based  on  discussion 
given  by  Sutton  (1953)  and  others  and  is  called  a "pragmatic"  approach. 
Here  6^  is  the  equivalent  potential  temperature  given  by 


0 = 
e 


and  A9  is  the  difference  between  the  values  of 
e 

the  lowest  GCM  level;  is  the  average  of  the  9 


9 at  the  surface  and 
e 

at  these  two  levels, 
e 
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The  main  purpose  of  using  Eq.  (lA)  is  to  be  able  to  establish 
(through  Q^)  a near  moist  adiabatic  lapse  rate  in  the  Im/est  km  as  is 
usually  observed  over  tropical  oceans,  and  at  the  same  time  to  allow  a 
steeper  lapse  to  exist  over  land  areas.  The  UK-I  GCM  uses 


A = 10^  m/s 


for  unstable  case,  i.e,, 


if  (A0  + e, ) > 0, 

el 


= 0 


for  other  cases. 


/ based  on  experience. 

^2  = 2°k) 


A limit  is  imposed  on  the  rate  of  heating  under  unstable  conditions 
such  as  when  extremely  cold  air  flows  over  warm  sea  such  that 


A 


A0  + £_ 

e 1 


1 


'a 


< 200  m/sec 


6 

e 


The  surface  humidity  mixing  ratio  is  calculated  on  the  assumption  that 
relative  humidity  at  the  ocean  surface  is  75  percent  and  that  at  the 
land  surface  is  50  percent.  This  enables  Eq.  (14)  to  provide  a real- 
istic discrimination  between  the  lapse  rates  over  a land  and  sea. 

Analogous  to  Eq.  (14),  the  latent  heat  flux  at  the  surface  is 
given  by 


E 

o 


(15) 


when  Cg,  A,  and  have  the  same  values  as  those  in  Eq.  (14). 

This  model  gives  unrealistic  results  since  strong  lapse  rates  may 
be  established  spuriously  without  instability.  As  a consequence, 
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one  may  have  dry  air  at  1 km  in  equilibrium  with  the  underlying  surface, 
which  IS  clearly  not  possible  over  desert  regions.  This  may  be  obviated 
by  distinguishing  between  different  types  of  land  surfaces. 

fietermination  of  Surface  Fluxes  in  UCLA-GCM  (3-Level  Version) 

The  parameterization  of  surface  eddy  fluxes  in  the  UCLA  GCM  is 
based  on  the  method  suggested  by  Deardorff  (1972).  This  method  deter- 
mines surface  fluxes  of  momentum,  heat,  and  moisture  as  functions  of 
stability  (bulk  Richardson  number) , surface  roughness  (Z  ) , and  the 

o ’ 

total  b,l,  depth  (h)  obtained  through  a rate  equation.  Deardorff  de- 
rived his  formulations  by  combining  the  explicit  formulations  for  the 
surface  layer  and  the  b.l.  deficit  formulations  for  the  layer  above 
the  surface  layer.  This  technique  enabled  him  to  eliminate  the  use  of 
the  values  of  the  parameters  at  the  so-called  anemometer  level  s,  and 

express  the  formulations  in  terms  of  the  bulk  properties  of  the  entire 

b.l. 

In  actual  practice,  the  UCLA  GCM  computes  the  surface  fluxes  in 
the  same  way  as  those  GCMs  which  use  Eq.  (12)  except  for  the  following 
differences : 


o The  surface  transfer  coefficients  C^^/C^  are  obtained  from 
Fxg.  6 which  shows  the  functional  dependence  of  these  coefficients  on 
the  bulk  Richardson  number  (Rig) , and  the  ratio  (h  - Z ) /Z  of  the  ef- 
fective b.l.  depth  to  the  roughness  parameter.  C^  is  assumed  to  be 
equal  to  C^.  The  figures  are  obtained  from  the  following  formulas  for 


(16) 


Here 


Ri. 


= _ 


(h  - z^)(e^  - 


C 8^. 

p M 


o The  model  uses  static  energy  instead  of  temperature  (0). 

o The  difference  A6  (Eq,  (12))  represents  in  this  case  the  ef- 

fective difference  between  the  vertically  averaged  static  energy  of  the 
b.l.  and  its  value  at  the  ground.  Similarly,  Aq  is  the  difference  be- 
tween the  value  of  q at  the  ground  and  the  mean  value  within  the  b.l. 

o This  method  also  considers  the  occurrence  of  saturated  b.l. 
through  its  entire  depth  as  well  as  the  presence  of  fog. 

o h,  the  top  of  the  b.l.,  is  not  fixed  but  instead  is  predicted 
by  a prognostic  equation. 

o Z = 0.45  m at  all  land  points,  and 
° -4 

= 2.3  X 10  m at  all  ocean  points. 


VERTICAL  PROFILES  OF  THE  TURBULENT  FLUXES 

As  indicated  earlier,  the  actual  method  of  parameterizing  b.l. 
in  GCM  depends  on  the  number  of  levels  (in  the  vertical)  incorporated 
within  the  b.l.  In  global  models  which  cannot  resolve  the  b.l.  ex- 
plicitly and  in  which  the  lowest  GCM  level  is  well  above  the  b.l.  top, 
only  the  surface  fluxes  are  of  any  significance.  In  some  global  models 
which  have  several  levels  within  the  b.l.,  the  vertical  structure  of 
the  b.l.  in  terms  of  the  vertical  profiles  of  turbulent  fluxes  can  be 
determined.  The  lowest  GCM  level  is  usually  placed  at  the  top  of  the 
surface  layer  (presumed  to  be  a constant  flux  layer)  and  the  other 
levels  are  placed  in  the  so-called  Ekman  layer  above.  We  have  described 
above  the  techniques  used  by  various  GCMs  to  evaluate  surface  fluxes. 

In  this  section  we  concentrate  on  the  determination  of  turbulent  fluxes 
within  the  Ekman  layer. 

The  formulations  generally  employed  to  compute  the  turbulent  fluxes 
in  the  layer  above  the  surface  layer  are 
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“>  dV 

stress  - T = pK^  — , 

heat  flux  . Q - -pCpKjj(||  - Y;,^)  . (17) 


moisture  = E = “P^s-  ^ > 

^ 32 


where  is  a counter-gradient  flux  (Deardorff,  1966)  and  K^,  and 

Kg  are  the  eddy  diffusion  coefficients  for  momentum,  heat,  and  moisture, 
respectively. 

The  actual  use  of  Eq.  (17)  in  GCMs  depends  on  the  representation 
of  the  b.l.  in  the  numerical  model,  and  these  can  be  considered  for  the 
following  three  cases : 

Case  I.  Implicit  b.l.  (NCAR,  GISS,  RAND  GCMs) 

NCAR,  GISS,  and  RAM)  GCMs,  which  do  not  incorporate  the  b.l.  ex- 
plicitly,  use  Eq.  (17)  to  compute  the  values  of  V^,  T^,  and  q^ , which 
are  required  to  evaluate  surface  fluxes  (Eq.  11).  The  method  of  com- 
puting  V^,  T^ , and  q^  is  based  on  the  assumption  that  there  is  no  flux 
divergence  in  the  b.l.  and  thus  the  surface  fluxes  can  be  equated  to 
the  fluxes  in  the  layer  from  level  s to  the  first  internal  GCM  level 
in  the  model.  For  example,  equating  surface  stress  to  the  stress  in 
the  Ekman  layer  we  get 


(18) 


Here  V^,  the  large-scale  wind  velocity  at  the  first  GCM  level,  is 
known  from  the  GCM,  and  if  C^and  arc  specified,  Eq.  (18)  can  be 
solved  for  the  unknown  V^.  Similar  sets  of  expressions  are  used  for 
obtaining  T^  and  q^,  which  are  then  used  to  determine  surface  fluxes 
for  heat  and  moisture  from  Eq,  (11). 
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Specification  of  K^/Kj^/K^ 

It  is  assumed  that  and  are  equal.  The  formulation  for 

is  given  by 


+ A^jl  - exp 


(~~y  ) 

!\3Z  XG/ 


for  neutral/unstable  case. 


h • 


36 

i.e. , ^ < T 


3Z  — CG  ’ 


(19) 


= + A 

1 + A,Ri  5 


for  stable  case,  i.e.,  — > 


where 


Ri  = Richardson  number 


f--Y  ) 
\3Z  XG/ 


^1’  ^2*  '^3*  ^5’  ^6  empirical  constants;  and  Ygg  is  the 

counter-gradient  heat  flux.  d is  a factor  introduced  to  reduce  the 
value  of  over  the  mountains  where  the  thickness  of  the  lowest  GCM 
layer  is  reduced  considerably. 

The  values  of  the  above  empirical  parameters  used  in  the  NCAR  GCM 
are  given  by 


10^  cm^/sec, 

10^  cm^/sec, 

1.2  X 10^  cm/°C, 
40, 

2 X 10^  cm^/sec. 
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5 X lo"^  °C/cm  at  Z = s , 


'CG 


1 X 10~^  °C/cm  otherwise  , 


d = 


1 - ^ for  H < AZ  , 


2 - “ for  AZ  < H < 2AZ  , 
AZ  _ — 


where  H is  the  height  of  orography.  The  GISS  GCM  uses  the  same  values 

for  A„,  A , A„ , A , A.,  A , and  A,  as  in  the  NCAR  GCM,  but  it  uses 
E 1 Z 3 H j D 

values  of  0 and  1 for  y__  and  d respectively.  The  RAND  GCM  uses 

LLr 


Case  II.  Explicit  b.l.  (GFDL/UK-II)  GCM 

Both  these  GCMs  place  the  lowest  GCM  level  at  the  top  of  the  sur- 
face layer.  GFDL  has  two  more  GCM  levels  up  to  the  top  of  the  b.l., 
and  the  turbulent  fluxes  in  the  layer  above  the  surface  layer  are  com- 
puted from  Eq.  (17)  through  finite  difference  formulations.  The  eddy 
diffusion  coefficients  are  assumed  to  be  the  same  for  momentum,  heat, 
and  moisture  fluxes  and  are  given  by 


3Z 


3Z 


(1  - a S)  for  unstable  case  ( < 0)  , 

c \ dZ 


(20) 


X / 9 G \ 

(1  + a S)  for  neutral/stable  case 


Here  £ is  the  mixing  length,  S is  a stability  parameter,  and  is  an 
empirical  constant  (a^  = 18).  The  b.l.  is  explicitly  included  in  the 
model  to  the  extent  that  a mixing  length  given  by  £ applies.  Table  2 
shows  the  values  of  £ and  S used  by  the  GFDL  and  the  UK-II  GCMs;  the 
formulations  have  also  been  used  by  Estoque  6i  Bhumralkar  (1969)  . 


Case  III.  Explicit  b.l.  Treatment:  UCLA  GCM  (3-Level) 

The  parameterization  technique  used  in  the  UCLA  GCM  is  essentially 
based  on  the  one  proposed  by  Deardorff  (1972)  and  is  designed  for  GCMs 


Table  2 


VALUES  FOR  £ AND  S IN  GCMS 


Model 


UK-II 


GFDL 


Total  Levels  Levels 


in  Vertical 


(i)  9(1969) 
(ii)18(1974) 


(ill)  9(1971)' 


in  b.l. 


1 

(850  m) 


30  meters 


(i)  30  meters 
(ii)  = 0.4Z  for  Z ^ s 


“ 0-4  ""g  _ for  s < Z < B 
s 


Z = 75  m ' 
s 

B = 2. 5 km, 


= 0 for  B < Z 


(iii) 


1 + 


0.4  Z 
1480  fZ 


0 (neutral) 


0 (neutral) 
0 (neutral) 


S = 


3Z 


3Z 


Remarks 


= 0 (Eq.  1/) 


Does  not  include 
vertical  heat  flux 


Also  Ygg  - 0 


Q- 

in  all  cases. 

B is  the  level  at 
which  £ = 0 


^ot  used  in  long-term  integrations,  but  used  by  Delsol  et  al.  (1971)  to  test  b.l.  parameterizations . 


oe- 
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that  do  not  resolve  the  b.l.  However,  the  method  involves  solving  pi^og- 
nostic  equations  for 


o the  height  of  the  b.l.  (h) , and 

o the  "strength”  of  the  discontinuity  in  momentum,  heat,  and 

moisture  at  the  b.l.  top  (AU,  AV,  A9,  Aq)  to  maintain  consis- 
tency between  the  b.l.  parameterization  and  the  GCM. 


The  first  item  has  been  discussed  in  considerable  detail  by  Bhumralkar 
(1974)  and  incorporates  the  entrainment  hypothesis  which  in  turn  de- 
pends, among  other  things,  on  the  flux  of  the  quantity  at  the  b.l.  top 
as  well  as  its  discontinuity  there.  Thus,  neglecting  radiative  effects. 


= - 


(21) 


where  is  the  heat  flux  at  the  top  of  the  b.l.,  and  A6  is  the  inver- 
sion "strength"  at  h. 

The  prognostic  equations  for  computing  the  "strength"  of  the  dis- 
continuity at  the  top  of  the  b.l.  require  evaluation  of  the  eddy  flux 
divergence  within  the  b.l.  (in  addition  to  advection  and  interaction 
with  cumulus  scale  processes) . Thus 


(22) 


It  can  be  seen  from  Eqs.  (21)  and  (22)  that  it  is  necessary  to  know 
the  vertical  distribution  of  eddy  fluxes  in  the  b.l.  This  is  also 
necessary  for  the  following  two  reasons: 


(a)  Since  the  b.l.  depth  is  variable  (and  determined  prognostic- 
ally),  the  b.l.  is  permitted  to  completely  fill  one  or  more 
GCM  layers  and  thus  the  eddy  fluxes  between  such  layers  must 
be  determined. 
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(b)  The  parameterization  method  considers  the  formation  of  new 

boundary  layers;  for  this  it  is  required  to  determine  the 

rate  of  conversion  of  the  buoyant  potential  energy  into  the 

turbulent  kinetic  energy  (B  = g/C  (Q/T)dZ),  which  in  turn 

p s 

depends  on  the  knowledge  of  vertical  distribution  of  the  flux 
of  heat  (say)  within  the  b.l. 

The  model  estimates  the  flux  of  heat  within  the  b.l.  from  the 
equation 


Q = q,  + [Q,  - V[^_ 


+ K E„  Sin 
N 


3Z 


(23) 


where  K is  a constant,  and  is  the  rate  of  entrainment  given  by  Eq. 
(21),  which  equals  0 at  the  surface  s and  the  b.l.  top  h.  Equation 
(23)  can  be  used  to  compute  eddy  flux  convergence  terms  at  h,  which 
are  then  used  to  determine  the  "strengths/discontinuities''  at  the  b.l. 
top  (reference  Eq.  22) . The  last  term  in  Eq.  (23) — the  mixing  term — is 
added  to  obviate  the  inconsistencies  which  arise  when  a simple  linear 
profile  of  h is  considered. 

An  equation  similar  to  (23)  can  be  written  for  moisture  flux  with- 
in the  b.l.  However,  at  this  time  (August /September  1974),  research 
is  continuing  and  no  conclusive  discussion  can  be  given. 

HORIZONTAL  DIFFUSION  TERMS  IN  GCMs 

There  is,  in  general,  some  disagreement  regarding  the  necessity 
to  include  explicit  horizontal  eddy  diffusion  terms  in  GCMs.  Until 
recently  there  was  a general  denial  of  its  physical  significance;  many 
felt  that  a parameterized  horizontal  diffusion  was  purely  mathematical 
gimmickry  to  suppress  computational  instability  (Smagorinsky,  1970) . 
However,  apart  from  the  development  of  differencing  methods  which  ap- 
parently obviate  the  need  for  these  terms,  the  experience  of  the  GFDL 
GCM  with  nonlinear  viscosity  designed  to  prevent  the  -5/3  power  law 

*For  details  see  Randall  and  Arakawa  (1974). 


near  the  limit  of  resolution  has  been  quite  satisfactory.  Whereas  the 
GFDL/NCAR/UK-I  GCMs  use  explicit  non-linear  lateral  diffusion  terms 
(Smagorinsky , 1963)  in  their  governing  equations,  the  RAND,  GISS,  and 
UCLA  GCMs  do  not.  However,  it  may  be  mentioned  that  a horizontal  aver- 
aging scheme  has  been  introduced  in  the  RAND  GCM  (and  also  the  UCLA  GCM) 
in  the  higher  latitude  regions.  The  UK-II  GCM  uses  a constant  diffusion 
coefficient  to  evaluate  lateral  diffusion. 


V.  EXPERIMENTS  TO  TEST  DIFFERENT  BOUNDARY-LAYER  FARAMETERIZATIQNS 
IN  AN  ATMOSPHERIC  GENERAL  CIRCULATION  MODEL 


NEED  FOR  EXPERIMENTS 

As  discussed  in  the  preceding  section,  several  approaches  to  the 
b.l.  parameterization  are  aimed  at  incorporating  the  b.l.  processes  in 
GCMs.  However,  in  the  absence  of  any  systematic  sensitivity  tests  it 
is  difficult  to  determine  the  specific  advantages  or  disadvantages  of 
each  of  the  approaches.  There  has  been  no  estimate  of  errors  of  param- 
eterization in  relation  to  the  basic  assumptions  of  the  theory  on  which 
a technique  is  based.  The  different  GCMs,  which  have  parameterized  the 
b.l.  processes  either  explicitly  or  implicitly,  do  not  specifically 
discuss  the  impact  of  the  particular  parameterization  technique  on  the 
global  results.  Thus,  there  has  been  no  attempt  to  test  the  degree  of 
acceptability  of  the  simplifications/assumptions  in  evaluation  of  tur- 
bulent fluxes  of  momentum,  heat,  and  moisture  in  GCMs.  In  fact,  there 
have  been  no  studies  even  to  determine  whether  the  incorporation  of 
b.l.  physics  in  a GCM  is  justified  or  not.  As  far  as  is  known  to  the 
author  at  this  time  only  Delsol  et  al.  (1971)  have  studied  (in  a numer- 
ical GCM)  the  effects  of  parameterization  of  various  processes  related 
to  the  b.l.  physics.  They  evaluated  the  comparative  effects  by  repeat- 
ing the  integrations  and  by  gradually  increasing  the  complexity  of  the 
model  in  terms  of  the  use  of 

o different  roughness  over  land  and  sea, 

0 Monin-Obukhov  treatment  of  the  surface  layer  turbulence, 
o stability  dependent  eddy  diffusion  coefficient  in  Ekman  layer, 
o diurnal  variation  of  insolation,  and 
o conduction  of  heat  into  the  soil. 

Their  results,  based  on  14  days  of  real-time  prediction,  showed  that 
the  sophistication  of  the  b.l.  physics  does  not  produce  a particularly 
large  effect  on  large-scale  prediction  until  about  7 days.  They  spec- 
ulate (without  evidence)  that  it  may  become  large  after  10  days. 


-35- 


Their  results  have  prompted  Pasquill  (1972)  to  suggest  that  (if  these 
are  true)  it  may  be  unneoessaTy  to  incorporate  much  of  the  detail  of 
the  b.l.  effects  in  GCMs. 

In  the  author’s  opinion,  this  conclusion  is  not  quite  justified 
because  the  experiments  of  Delsol  et  al.  (1971)  are  not  conclusive  due 
to  various  factors.  For  example: 

(1)  They  assume  a mixing  length  of  30  m,  and  yet  their  first  GCM 
level  is  at  75  m,  the  next  at  640  m,  and  so  on.  This  does 
not  enable  the  model  to  resolve  the  vertical  eddy  size  in 
both  the  constant  flux  layer  and  the  Ekman  layer. 

(2)  The  formulations  used  in  the  Ekman  layer  are  not  quite  valid 
in  low-latitude/equatorial  regions. 

(3)  No  comparison  with  observed  data  was  attempted. 

In  view  of  the  above,  and  because  of  the  increasing  importance  of 
the  use  of  GCMs  for  climatic  studies,  it  is  evident  that  a comparison 
of  various  b.l.  parameterization  schemes  is  desirable  in  order  to  se- 
lect the  "best”  approach.  This  view  has  also  been  expressed  at  various 
international  meetings  held  to  discuss  the  problem  of  b.l.  parameter- 
ization (GARP  Publications  Series,  No.  8,  October  1972,  and  Report  No. 
5,  July  1974,  of  the  GARP  Programme  on  Numerical  Experimentation). 

PARAMETERS  INVOLVED  IN  COMPARISON  TESTS 

Before  any  experiments  are  performed  to  compare  different  b.l. 
parameterization  schemes  in  a GCM,  it  is  necessary  to  specify  param- 
eters or  groups  of  parameters  in  relation  to  which  the  comparisons 
should  be  performed.  It  appears  that  it  may  be  useful  to  integrate  a 
GCM  by  including  different  b.l.  parameterization  schemes  (one  at  a 
time)  and  then  compare  these  specified  parameters  as  functions  of  the 
parameterization  schemes.  For  a consistent  and  meaningful  comparison 
it  is  necessary  that  for  each  integration  the  general  prognostic  scheme 
and  initial  conditions  should  remain  fixed.  The  comparison  may  be  per- 
formed, over  the  globe,  with  respect  to  the  following  parameters  under 
all  stability  conditions: 
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(1)  Values  of  the  turbulent  surface  fluxes  of  moinentum,  heat,  and 
moisture. 

(2)  Vertical  velocities  at  the  top  of  the  b.l. 

(3)  Height  of  the  b.l. 

(4)  Global  heat  flux  distribution  in  the  vertical. 

(5)  Vertical  profiles  of  the  turbulent  fluxes  as  well  as  of  the 
meteorological  variables  within  the  b.l. 

(6)  Variability  of  surface  parameters  (roughness,  soil  moisture 
content,  etc.)  at  subgrid  scale.  This  involves  use  of  com- 
monly used  aerodynamic  roughness  length  vs  the  effective 
roughness  (Fiedler  and  Panofsky,  1972). 

It  is  envisaged  that  in  order  to  determine  the  applicability  of 
a b.l.  parameterization  scheme  to  a GCM,  the  results  pertaining  to  the 
above  quantities  will  be  compared  between  the  different  models  them- 
selves, as  well  as  with  the  observational  data  for  some  situations  for 
which  there  are  sufficient  data  available. 

Another  important  factor  which  could  influence  the  determination 
of  the  quality  of  a parameterization  scheme  (in  a GCM)  is  the  real  time 
du'pat'ton  of  i.ntegTat'ion  of  the  GCM.  Zilitinkevich  (1972)  has  suggested 
that  it  may  be  reasonable  to  compare  results  of  different  integrations 
for  3 to  4 days.  However,  as  stated  earlier,  the  experience  of  Delsol 
et  al.  (1971)  indicates  that  even  7 days  of  integration  may  not  be  suf- 
ficient to  bring  out  significant  differences.  Thus  it  appears  that  the 
integrations  may  have  to  be  performed  for  at  least  a two-week  period 
or  even  longer.  It  may,  however,  be  noted  that  due  to  lack  of  exper- 
ience, either  the  specification  of  integration  period  or  the  list  of 
"comparable"  parameters  (above)  is  not  final  and  these  may  have  to  be 
modified . 

PROPOSED  EXPERIMENTS 

We  propose  to  perform  a series  of  experiments  to  compare  some  b.l. 
parameterization  schemes  on  the  basis  of  discussions  given  in  the  pre— 
two  subsections . We  expect  to  use  an  existing  version  of  the 
Rand  2-level  GCM  as  a tool  to  perform  these  experiments.  It  is  evident 
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that,  due  to  the  coarse  vertical  resolution  of  the  Rand  GCM,  we  will 
not  be  able  to  determine  the  vertical  structure  of  the  b.l.  in  respect 
to  the  profiles  of  the  fluxes  as  well  as  of  the  meteorological  vari- 
ables. Consequently,  the  emphasis  will  have  to  be  on  the  determination 
and  comparison  of  surface  fluxes  with  observations. 

The  proposed  experiments  may  involve  comparison  of  parameterization 
schemes  which  are  based  on  similarity  theory  as  well  as  K-theory  con- 
cepts. For  a given  configuration  of  the  GCM  prognostic  scheme  and  in- 
itial conditions,  we  propose  to  integrate  the  model  by  incorporating 
the  following  succession  of  parameterization  schemes: 

(1)  K-theory  scheme 

The  formulations  similar  to  the  one  used  by  the  NCAR/GISS 
GCM  may  be  used.  This  may  be  performed  for  different  values  of 
the  eddy  coefficient  K and  of  the  drag  coefficient. 

(2)  Similarity  theory  scheme 

(a)  We  can  use  the  formulations  in  their  original  form  which 
consider  the  scale  height  of  b.l.  to  be  given  by  U^/f.  The  dis- 
advantage of  using  f in  low  latitude/equatorial  regions  can  be  ob- 
viated by  specifying  a critical  latitude  and  using  f for  this 
across  the  equatorial  region  between  the  critical  latitudes  on 
either  side  of  the  equator. 

(b)  Instead  of  the  height,  h,  of  the  b.l.  can  be  com- 

puted by  a prognostic  equation.  The  latter  involves  use  of  en- 
trainment hypotheses  and  can  be  patterned  after  Deardorff  (1972). 

(c)  Instead  of  U^/f  and  a prognostic  equation  for  h which 
may  give  rise  to  singularities  under  certain  conditions,  we  can 
also  test  the  use  of  "interpolation"  prognostic  equations,  such 
as  that  proposed  by  Deardorff  (1974) , which  has  been  designed  to 
overcome  the  singularities. 

(d)  The  above  three  procedures  involve  a given  set  of  so- 
called  "universal"  functions.  It  may  be  found  necessary  to  exper- 
iment with  other  sets  of  values  for  these  functions. 

(e)  The  effect  of  the  use  of  the  roughness  parameter  (z^) 
can  also  be  investigated.  In  this  work  it  may  be  desirable  to 
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test  the  formulations  suggested  by  Zilitinkevich  (1969)  and  Clark 
(1970).  The  concept  of  effective  roughness  (Fiedler  and  Panofsky, 
1972)  can  also  be  investigated. 

(3)  Variable  boundary- layer  thickness 

In  addition  to  the  above  we  can  perform  an  experiment  using 
the  parameterization  scheme  developed  by  UCLA  for  their  3-level 
GCM.  In  the  Rand  2-level  GCM,  however,  we  can  treat  only  that 
case  which  considers  an  entirely  sub-grid  scale  boundary  layer. 
This  is  because  in  this  model  it  may  be  more  realistic  to  assume 
that  the  b.l.  depth  is  always  confined  to  the  lower  layer  (sur- 
face to  600  mb)  and  is  thus  not  resolved  by  the  GCM, 

The  results  and  conclusions  of  such  comparisons , together  with 
specific  details  in  respect  of  the  parameterization  schemes  used,  will 
form  the  subject  matter  of  a subsequent  report. 
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VI.  CONCLUDING  REMARKS 


It  is  recognized  that  we  can  hope  to  make  a quantitative  assessment 
of  the  general  circulation  and  climatic  changes  only  through  the  use  of 
reasonably  calibrated  numerical  models.  With  this  end  in  view,  there 
have  been  concerted  efforts  to  improve  the  fidelity  of  GCMs  by  improv- 
ing (among  other  physical  processes)  the  parameterization  of  boundary- 
layer  (b.l.)  processes,  by  using  schemes  that  are  not  unique  and  can  be 
approached  in  a variety  of  ways. 

At  present  there  are  several  approaches  to  the  b.l.  parameteriza- 
tion in  GCMs , and  it  is  fair  to  expect  that  the  number  of  these  formu- 
lations will  increase  in  time.  However,  so  far  there  have  been  no 
studies  aimed  at  the  comparison  of  the  various  parameterization  schemes 
in  GCMs  in  order  to  select  the  best  approach  and  to  estimate  errors  of 
the  parameterizations.  It  is  evident  that  the  limitations  associated 
with  the  theories  used  in  parameterizations  will  have  a controlling  in- 
fluence on  the  results.  For  example,  in  case  of  similarity  theory,  it 
is  not  clear  to  what  extent  the  assumptions  of  this  theory  as  regards 
horizontal  homogeneity  of  the  underlying  surface  (in  the  steady  state) 
must  be  satisfied  to  permit  utilization  of  the  conclusions  of  the  the- 
ory. Also,  the  heterogeneity  of  the  underlying  surface  is  currently 
taken  into  account  by  using  the  roughness  parameter  (Z^)  as  a function 
of  the  character  of  the  underlying  surface.  It  is  yet  to  be  determined, 
however,  whether  this  is  the  only  possible  kind  of  heterogeneity  on  a 
sub-grid  scale.  The  use  of  a rate  equation  to  obtain  the  height  of 
the  b.l.  has  also  to  be  investigated  before  the  method  is  universally 
applied  in  practice. 

What  we  have  suggested  in  this  review  is  a systematic  approach  to 
determine  the  most  suitable  b.l.  parameterization  for  the  Rand  GCM. 

Through  comparison  of  the  simulation  responses  of  the  GCM,  the 
experiments  can  be  designed  to  determine  the  relative  role  and  the 
sensitivity  of  the  various  formulations. 
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